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Abstract — Floods result from excessive runoff quantities 
and rates, which usually cause costly damage to properties and 
sometimes loss of lives. Flood routing is the spatio-temporal 
monitoring of flood wave as it traverses an area. Flood routing 
techniques are generally classified as hydrologic and hydraulic. 
Flood routing of Asa River was conducted using kinematic 
wave method, a hydraulic technique. The Asa River is an 
important river in Ilorin, the capital of Kwara State, Nigeria. 
The United States Soil Conservation Service (SCS) method was 
used in developing the synthetic unit hydrograph using the 
computed peak discharge of 2.44 m 3 /s and time to peak of 0.213 
hr. Kirpich’s equation was used to compute the time of 
concentration of the watershed as 0.32 hr. 

Index Terms — Flood routing, Kinematic wave, Hydrograph 

I. INTRODUCTION 

Flooding may result from the volume of the water within a 
body of water, such as a river or lake, exceeding the total 
capacity of its bounds, with the result that some of the water 
flows outside the catchment. It can also occur in rivers, when 
the strength of the river is so high it flows out of the river 
channel, particularly at bends or meanders (Khatibi, 2012). 

The sudden change in the inflow due to occurrence of heavy 
rain in the upstream usually leads to the unsteadiness in a 
river as well as breaches in the embankment system. This 
may results to the collapse of a river valley structure 
(Sudarshan and Arup, 2013) 

Framji and Garg (1976) defined flood as a relatively high 
flow in a river above normal, also the inundation of low land 
which may result from the high flow. Fread (1981) and 
Linsley et al., (1989) defined flood routing as a mathematical 
method for predicting the changing magnitude and celerity 
of a flood wave as it propagates through rivers or reservoirs. 
Generally, flood routing methods are categorized into two 
broad, but somewhat related applications, namely reservoir 
routing and open channel routing (Lawler, 1964). These 
methods are frequently used to estimate inflow or outflow 
hydrographs and peak flow rates in reservoirs, river reaches, 
farm ponds, tanks, swamps and lakes (SCS, 1972; Viessman 
et al., 1989; Smithers and Caldecott, 1995). 
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Flood routing is an important technique used in the design of 
flood protection measures in order to estimate how the 
proposed measures will affect the behaviour of flood waves in 
rivers so that adequate protection and economic solutions can 
be found (Wilson, 1990). In practical applications, two steps 
are involved in the prediction and assessment of flood level 
inundation. To estimate the flow-gauging station to a 
downstream location, a flood routing model is employed. 
Then the flood hydrograph is input to a hydraulic model in 
order to estimate the flood levels at the downstream site 
(Blackburn and Hicks, 2001). Routing is used to predict the 
temporal and spatial variations of a flood wave as it traverses 
a river reach or reservoir, or it can be employed to predict the 
outflow hydrograph from a watershed subjected to a known 
amount of precipitation (Viessman et al., 1989). 

Flood-routing procedures may be classified as either 
hydrological or hydraulic (Choudhury et al., 2002). 
Hydrological methods use the principle of continuity and a 
relationship between discharge and the temporary storage of 
excess volumes of water during the flood period (Shaw, 
1994) while hydraulic methods of routing on the other hand 
use both the equation of continuity and equation of motion, 
customarily the momentum equation (Michael and Ojha, 
2006). This particular form utilizes the partial differential 
equations for unsteady flow in open channels. It adequately 
describes the dynamics of flow more than does the hydrologic 
routing techniques. 

Hjelmfelt (1984) stated that kinematic wave equation can be 
expressed in the form of a convolution integral. The kernel of 
this convolution integral can then be compared to a unit 
hydrograph. Kinematic wave equations are non-linear and 
hyperbolic; their exact solution is found using the method of 
characteristics (Wang and Tartakovsky, 2012). 

Field (1982) proposed using the method of characteristics to 
solve the kinematic wave equation to produce the discharge 
hydrograph. In his work, he demonstrated that the model is a 
generalization of several other popular models and typical 
values of some catchment parameters were suggested for 
their application to design. 

Numerical technique for solving the kinematic wave 
equation requires the initial values for the entire domain of 
the distance travelled by the flood as well as the upstream 
boundary conditions for the time of travelled. At each of the 
grid point, finite difference method can be adopted to obtain 
the solutions of the kinematic wave equation (Hossain et al., 
2013). The development of the kinematic wave flow routing 
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equation and certain of the fundamental conditions necessary 
for its successful application are discussed by Smith (1980). 

The objectives of this research work were to route an inflow 
hydrograph through a given river reach and obtain its 
outflow hydrograph and also to determine the amount of 
attenuation in hydrographs between the inflow and outflow at 
the end of a routing reach using kinematic wave method. 


II. METHODOLOGY 

A. The Study Area 

Asa River has its source from Oyo State, Nigeria and it flows 
through Ilorin, Kwara State, Nigeria in a South-North 
direction forming a dividing boundary between Eastern and 
Western Ilorin. It is about 56 km long, with a maximum 
width of about 100 m (at the dam site). Asa River has its 
estuary at River Awon, which is one of the tributaries of 
River Niger, at 12,200 m North of Ilorin. It is joined by River 
Oyun to the East and to the West by River Imoru. Afidikodi, 
Ekoro, Obe are among the earliest tributaries of Asa River 
while its tributaries in Ilorin include River Agba, Aluko, 
Atikeke, Mitile, Odota, Okun, and Osere. (Ojo, 1998; 
Ibrahim et al., 2013). 

The Asa catchment in Ilorin, North Central Nigeria is 
located between latitudes 8°36'N and 8°24'N and longitudes 
4°36'E and 4°10'E. Its total catchment is 1037 km 2 at the 
confluence and lies within Kwara State and Oyo State, in 
Nigeria with about one third the basin area in Oyo State 
(Okekunle, 2000). Asa River is a very significant source of 
water in terms of economic, agricultural and environmental 
purposes in the city as it is used in homes and industries 
(Ahaneku and Animashaun, 2013). There are farmlands, 
residential and industrial buildings along the bank of the 
river downstream. The catchment is formed by a ridge of 
hills rising to almost 580 m above sea level and in most 
places the catchment is a gently undulating plane. The soils 
in the catchment area are as a result of weathering of parent 
rocks, which support vegetation along the source of the river 
(Salawu, 1987). 

The Asa Dam is located at a point about 5 km south of Ilorin 
across Asa River. It is a composite dam of earth embankment 
at its extreme ends, with a central spillway followed to the 
right by a mass concrete non-overflow gravity section. The 
dam is 597 m long and 27 m high at its deepest section and a 
crest width of 6 m. A stilling basin of the entire width of the 
spillway dissipates the energy of the spill flow thereby 
preventing erosion of the stream bed (Kwara State Water 
Corporation, 1985). 

B. Development of Synthetic Unit Hydro graph 

The SCS, Soil Conservation Service - Natural Resources 
Conservation Service, NRCS method was used for 
computing unit hydrograph needed to develop the storm 
hydrograph. The time to peak and the peak discharge were 
estimated using the area of the watershed, slope and length of 
travel. 


Watershed area, A = 2.5 km 2 

Length of travel, L = 0.289 km 

Watershed slope, S = 0.0014 (Olaniran, 1983) 


The time of concentration, t Ci was computed using Kirpich’s 
equation 

L 

t c = 0.066[— ] 077 (1) 

\S 


where: 


L = length of travel (km) 

S = watershed slope m/m 
t c = 0.32 hr 

Duration of excess rainfall, D was calculated using equation 

( 2 ) 

D = 0.133 t c (2) 


D = 0.042 hr 

Time to peak, t p was then calculated using equation (3) 

(*c + D) 

(3) 


t P = 

1,7 

tp= 0.213 hr 


Peak discharge was calculated from equation (4) 

(0.20SAQ) 

q P = 7 (4) 


q p = 2.44 m 3 /s 


The calculated values for peak discharge, q p , and time to 
peak, t p , were used to obtain points for unit hydrograph flow 
rate and its corresponding time. 


Table 1: Hie sectional geometry and other hydraulic properties of Asa River 


Distance from Width, B 
Initial Point (m) (m) 

Hydraulic 
Depth. Y (m) 

Area, A 
(m^) 

Wetted Perimeter 
(m) 

Hydraulic 
Radius, R(m) 

Discharge 

(m J /s) 

1.5 

1.75 

0.73 

1.2S 

3.2! 

0.52 

0.73 

3.5 

2.25 

0.75 

1.69 

3.75 

0.56 

0.63 

6.0 

2.00 

0.76 

1.71 

3.52 

0.62 

0.51 

SO 

2.00 

1.2S 

2.56 

4.56 

0.7S 

0.95 

10.0 

2.00 

1.55 

3.10 

5.10 

0.S7 

1.57 

12.0 

2.00 

2.50 

5.00 

7.00 

1.11 

3.08 

H.O 

2.00 

2.56 

5.12 

7.12 

1.12 

3.97 

ISO 

2.00 

2.13 

4.26 

6.26 

1.03 

3.32 

20.0 

2.00 

2.S0 

5.60 

7.60 

1.17 

3.11 

22.0 

2.00 

2.S0 

5.60 

7.60 

1.17 

3.14 

240 

2.00 

2.26 

4.52 

6.52 

1.06 

2.65 

26.0 

2.00 

1.49 

2.9S 

4.9S 

0.S5 

1.92 

2S.0 

2.00 

1.34 

2.6S 

4.6S 

O.SO 

1.93 

30.0 

2.00 

1.55 

3.10 

5.10 

0.S7 

2.16 

32.0 

2.00 

1.34 

2.6S 

4.6S 

0.80 

1.74 

34.0 

2.00 

1.43 

2.S6 

4.S6 

0.S3 

2.00 

36.0 

1.9 

1.22 

2.32 

4.34 

0.74 

1.54 


Source: Kwara State Water Corporation, Ilorin, Kwara State. 


C. Development of Storm Hydro graph 

Storm hydrographs of 10-yr, 24-hr; 25-yr, 24-hr; 50-yr, 
24-hr; and 100-yr, 24-hr were developed from the synthetic 
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unit hydrograph using convolution procedures (Wanielista, 
1990), with rainfall data from Ogunlela (2001). The 
cumulative runoff and the watershed potential storage were 
estimated from: 


Q= (P - 0,25) 2 

for P > 0.2S 



0.8S) 



where: 


(Q = 0 for P < 0.2S) 
(25 400) 

S = [~ -] - 254 

CN 



Q = cumulative runoff (mm) 

S = potential storage (mm) 

P = cumulative precipitation (mm) 

CN = curve number 

A curve number of 71 was obtained for the Asa river 
watershed from its hydrologic, soil and land - use 
characteristics. 


Holding and Stephenson (1995) stated that for numerical 

Ax 

stability, the grid ratio - — must be greater than or equal to 

At 

the wave celerity, c, where Ax = distance interval and At = 
time step. 

Kinematic wave velocity, U was estimated using (Hussain 
and Ferdous, 2013): 

u = v + 7g^ (7) 

where: 

V = average velocity of flow; 0.588 m/s (from Table 

1) 

g = acceleration due to gravity; 9.81 m/s 2 
y = average depth of flow; 1.701 m (from Table 1) 
U = 4.67 m/s 


Time required for a small kinematic wave to traverse a 
distance interval Ax, K, was calculated using; 

Ax 

K = (8) 

u 


where: 

Ax = change in length of the river (at any 
interval). (100 m was used) 

U = kinematic velocity; 4.67 m/s as previously 

calculated. 

K = 21.4 secs ~ 21 secs 

The area at any given point interval Ax was calculated using; 
At K 


A2,2 — 


I 


K+ At 


A2,2 + 


K+ At 


Ai,2 + 


(9) 


U (K + At) 
where: 

A 12 is the original area at point where Ax = 0; Ai > 2 = 
61.314 m 2 

A 2 ,i is the assumed area between A u and A 2 , 2 ; A 2 ,i = A 2 ,2 
- 0.02 


A 2 ,2 is the new area at every length interval, Ax 
At is the time step which equals to (K - t) and At < Ax/U 
(Courant condition). 

And t = 2 secs. Therefore, At = (K - 2) in secs. 

The outflow was computed using (Barati et al., 2012): 


Q2,2 - 

Ax 
+ 9At 


-Ax 

9At 

A2,l 


A2,2 + [ 



0 

_ (Ql,l~ Q2,l) + Ql,2 ] 

( 10 ) 


where: 


Qi i = the first inflow 

Qi 2 = second inflow 

Q 2 ,i = first outflow 

Q 2 ,2 = second outflow 

0 = weighting factor (0 = 0.75) 

Since A 2 ,i(assumed area) is 0.02 greater than Ai j2 (area at 
point where Ax = 0), therefore A 2 ,2 (new area at Ax interval) 
is 0.02 greater than A 2 ,i(assumed area). 

Before routing, it was assumed that the initial inflow equaled 
to the initial outflow. After the outflow Q 2 ,2 was calculated, it 
was used as I in equation (9) to get new A 2 , 2 - The computed 
A 2 ,2 was applied in equation (10) and subsequently used as 
Ai j2 in equation (9) for the next iteration to get new A 2 ,2 
which was used in equation (10) to get new Q 2 ,2 and so on. 


III. RESULTS AND DISCUSSION 

A. Results 

The results obtained from the whole procedures and all the 
parameters involved in routing the gauged Asa River are 
summarized in Table 2. The synthetic unit hydrograph peak 
discharges ranged from 0.01m 3 /s to 2.44m 3 /s while the time 
to peak ranged from 0.02 hr to 1.07 hrs. 

The 10-yr, 24-hr, 25 -yr, 24-hr, 50-yr, 24-hr and 100-yr, 
24-hr storm hydrographs developed gave peak flows of 73.47 
m 3 /s, 97.01 m 3 /s, 113.41 m 3 /s and 134.61 m 3 /s, respectively 
with time to peak of 0.213 hr. 

The Inflow and Outflow Hydrographs; Synthetic Unit 
Hydrograph; 10 yr, 24 hr Design Storm Hydrograph; 25 yr, 
24 hr Design Storm Hydrograph; 50 yr, 24 hr Design Storm 
Hydrograph and 100 yr, 24 hr Design Storm Hydrograph for 
Asa River are shown in figures 1,2, 3, 4, 5 and 6 respectively. 

B. Discussion 

The time of concentration obtained from the estimated 
watershed parameters of the catchment 0.32 hr as shown in 
Table 2, was used to compute the time peak and peak 
discharge requires developing a synthetic unit hydrograph 
using the dimensionless hydrograph ratio developed by U.S. 
Department of Agriculture Soil Conservation Service (SCS) 
(1972), now the Natural Resources Conservation Services 
(NRSC). The corresponding values of 2.44 m 3 /s and 0.213 hr 
were obtained for peak discharge and time to peak 
respectively. The developed SCS unit hydrograph was used 
to construct storm hydrographs of 10 yr, 24 hr; 25 yr, 24 hr; 
50 yr, 24 hr; 100 yr, 24 hr. 
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The peak discharge before routing was 3.97 m 3 /s at 1.8 hr 
and the peak discharge after routing was 3.87 m 3 /s at 1.8 hr 
indicating 0.1 m 3 /s attenuation. 

Also, the results show that there was some amount of water 
stored within the reach as well as a delay in the flow of water 
within the reach which indicates the occurrence of 
attenuation. 

The effect of the discharge will be felt at the upstream of the 
catchment towards the downstream. There were losses due to 
seepage, evaporation, inflow into minor tributaries and 
underground reservoirs. These losses affected the attenuation 
and storage within the reach. 

IV. CONCLUSIONS 

The routing of Asa River was carried out using Kinematic 
wave method. This method exhibited some level of accuracy 
as shown from the results in Table 2 when compared to some 
other methods of routing. Kinematic wave is one of the 
hydraulic methods of routing that takes into consideration 
watershed characterics and other hydraulic properties of the 
channel. Moreover, this method of flood routing is tedious 
and time consuming especially when all the required data are 
not available. 


V. Recommendations 

For proper and effective control of flood in the Ilorin 
metropolis, there should be proper maintenance of all 
hydraulic structures along the river, and due to the rapid 
development in the town there should be no erection of any 
structures along the flood plains. There should be proper 
channelization and trenching of the river from adequate 
monitoring of the river to forecast flood. Refuse disposal into 
the river should be prohibited and be backed up with 
legislation. 


Table 2: Summary of Results 


Parameters 

Results 

Watershed (Estimated) Area, A 

2.5 km 2 

Length of travel, L 

0.289 km 

Watershed slope, S 

0.0014 

Curve Number 

71 

Inflow Svnthetic unit hvdrograph 


Time of concentration, T c 

0.32 hr 

Duration of Excess rainfall, D 

0.042 hr 

Time to peak, T p 

0.213 hr 

Peak discharge, q p 

2.44 m 3 /s 

Kinematic Parameters 


Kinematic wave velocity, U 

4.67 m/s 

Time step, At 

19 sec 

Time required to traverse, K 

21 sec 

Weighting factor, 0 

0.75 

Change in length, Ax 

100 m 



Figure 1 : Inflow and Outflow Hydrographs for Asa River 



Figure 2: Synthetic Unit Hydrograph for Asa River 



Figure 3: 10 yr, 24 hr Design Storm Hydrograph 



Figure 4: 25 yr, 24 hr Design Storm Hydrograph 


16 


www.erpublication.org 


International Journal of Engineering and Technical Research (IJETR) 

ISSN: 2321-0869, Volume-2, Issue-3, March 2014 



Figure 5: 50 yr, 24 hr Design Storm Hydrograph 



Figure 6: 100 yr, 24 hr Design Storm Hydrograph 


REFERENCES 

[1] Ahaneku, I. E. and Animashaun, I. M. (2013). Determination of water 
quality index of river Asa, Ilorin, Nigeria. Advances in Applied Science 
Research. Vol. 4(6): 277-284. 

[2] Blackburn, J. and Hicks, F. E. (2001). Combined Flood and Flood level 
Forecasting. Proceeding 15 th Can. HydroTech. Conference. Victoria, 
British Columbia, Canada. 345-350. 

[3] Banti, R., Rahimi, S. and Akbari, H. G. (2012). Analysis of dynamic 
wave model for flood routing in natural rivers. Journal of Water Science 
and Engineering. Vol. 5(3): 243-258. 

[4] Choudhury, P., Shrivastava, R. K. and Narulkar, S. M. (2002). Flood 
routing in river networks using equivalent Muskingum inflow. Journal 
of Hydrology, 7 (6) 413-419. 

[5] Engineering Group (2004). The Ohio Department of Natural Resources, 
Division of Water, Dam safety: Probable Maximum Flood, Columbus. 

[6] Field, W. G. (1982). Kinematic Wave Theory of Catchment Response 
with Storage. Journal of Hydrology, Vol. 55: 279 - 301. 

[7] Framji, K. K. and Garg, B. C. (1976). Flood Control in the World. 
McGraw -Hill Book Company Inc. 

[8] Fread, D. F. (1981). Flood routing: A synopsis of past, present, and 
future capability. Proceeding International Symposium Rainfall-Runoff 
Modelling. Mississippi, USA. 521-541. 

[9] Hjelmfelt, Jr, A. T. (1984). Convolution and the Kinematic Wave 
Equations. Journal of Hydrology, Vol. 75: 301 - 309. 

[10] Holding, A. P. and Stephenson, D. (1995). Finite - Difference 
Formulations of Kinematic Equations. Journal of Hydraulic 
Engineering, Vol. 121(5):423 - 426. 

[lljHossain, M. M. and Ferdous, J. (2013). Solution of kinematic wave 
equation using finite difference method and finite element method. 
Global Journal of Science Frontier Research Mathematics and Decision 
Sciences. Vol. 13(6): 25-35. 

[12]Hossain, M. M., Ferdous, J. and Hossain, M. M. (2013). Sensitivity 
Analysis for the Physiographic Parameters of Overland as well as 
Channel Flow Elements Based on Kinematic Wave (Kw) Theory. 
Global Journal of Sciences Frontier Research Mathematics and Decision 
Sciences, Vol. 13(17): 29 -42. 


[13] Ibrahim, K. O., Okunlola, I. A. and Abdurrahman, A. (2013). Trace 
metal indices in the characterization of hydrogeochemical condition of 
surface water along Asa River, Ilorin, Kwara State, Nigeria. 
International Journal of Geology, Earth and Environmental Sciences. 
Vol. 3(l):29-35. 

[14] Khatibi, B. M. (2012). The study and estimate of floods in meandering 
river channels according to changes in the channel bending (case study 
Aji-Chay river in the North West of Iran). Geography and 
Environmental Hazards. Vol. 3:5-6. 

[15] Kirpich, Z. P. (1940). Time of concentration of small Agricultural 
watersheds. Civil Engineering (N. Y.), Vol. 10(6): 362. 

[16] Kwara State Water Corporation, (1985). The Sectional Geometry and 
other Hydraulic Properties of Asa River, Ilorin, Kwara State. 

[17] Fawler, E. A. (1964). Flood routing. US Army Corps of Engineers (ed.) 
Handbook of Applied Hydrology, Section 25-11 River Division. US 
Army Corps of Engineers, Ohio, USA. 35-58. 

[18] Finsley, R. K., Kohler, M. A. andPaulhus, J. F. H. (1989). Hydrology for 
Engineers. McGraw-Hill New York, USA. 

[19] Michael, A. M. and Ojha, T. P. (2006). Principles of Agricultural 
Engineering, Volume II. M/S Jain Brothers, New Delhi. 

[20] Ogunlela, A. O. (2001). Predicting Effects of Fand Use Changes on 
Runoff using the Curve Number Method. Nig. Journal of Tropical 
Engineering Vol. 2(1) 25 - 32. 

[21] Ogunlela, A. O. and Kasali, M. Y. (2002). Evaluation of Four Methods 
of Storm Hydrograph Development for an Ungauged Watershed. 
Nigerian Journal of Technological Development, Vol. 2 (1) 15 - 24. 

[22] Ojo, O. D. (1998). Sediment Transport along Kunrun Stream and Asa 
River, Ilorin. B.Eng. Project Report, Department of Agricultural 
Engineering, University of Ilorin, Ilorin, Nigeria. 

[23] Okekunle, R. C. (2000). Flood Routing along Asa River (Ilorin) using 
Kinematic Wave and Convex Routing Methods. B.Eng. Project Report, 
Department of Agricultural Engineering, University of Ilorin, Ilorin, 
Nigeria. 

[24] Olaniran, O. J. (1983). Flood Generating Mechanism at Ilorin, Nigeria. 
Journal of Geography Vol. 7(3): 27 1 - 277. 

[25] Salawu, K. F. (1987). Flow Augmentation for Improving Assimilation 
Capacity of Asa River, M. Eng. Thesis. Department of Civil 
Engineering, University of Ilorin, Ilorin, Nigeria. 

[26] Shaw, E. M. (1994). Hydrology in Practice. T J Press (Padstow) FTD, 
Cornwall, UK. 

[27] Smith, A. (1980). A Generalized Approach to Kinematic Flood Routing, 
Journal of Hydrology, vol. 45: 71 - 89. 

[28] Smithers, J. C. and Caldecott, R. E. (1995). Hydrograph routing. In: 
Schulze R.E. (ed.) Hydrology and Agrohydrology: A Text to 
Accompany the ACRU 3.00 Agrohydrological Modelling System. WRC 
Repot No TT 69/95. Water Research Commission, Pretoria, RSA. 
AT13-1 - AT13-16. 

[29] SCS. Soil Conservation Service, (1972). Flood Routing. Soil 
Conservation Service, US Department of Agriculture, Beltsville, USA. 

[30] Sudarshan, P. and Arup, K. S. (2013). Hydrodynamic flood routing 
considering piedmont zone. International Journal of Civil and Structural 
Engineering. Vol. 3(3): 464-474. 

[31] Viessman, W. Fewis, G. F. Knapp, J. W. and Harbough, T. E. (1989). 
Introduction to Hydrology. Harper and Row Publishers, New York, 
USA. 

[32] Wang, P. and Tartakovsky, D. M. (2012). Uncertainty quantification in 
kinematic-wave models. Journal of computational Physics. Vol. 
231:7868-7880. 

[33] Wanielista, M. P. (1990). Hydrology and Water Quality Control. John 
Wiley and Sons. Inc. USA. 

[34] Wilson, E. M. (1990). Engineering Hydrology. MacMillan, Hong Kong, 
China. 


17 


www.erpublication.org 


